Standard Quantum Limit for Probing Mechanical Energy Quantization 
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We derive a standard quantum limit for probing mechanical energy quantization in a class of 
systems with mechanical modes parametrically coupled to external degrees of freedom. To resolve 
a single mechanical quantum, it requires a strong-coupling regime — the decay rate of external 
degrees of freedom is smaller than the parametric coupling rate. In the case for cavity-assisted 
optomechanical systems, e.g. the one proposed by Thompson et al. [ll], zero-point motion of the 
mechanical oscillator needs to be comparable to linear dynamical range of the optical system which 
is characterized by the optical wavelength divided by the cavity finesse. 

PACS numbers: 



Introduction. — Recently, significant cooling of mechan- 
ical modes of harmonic oscillators has been achieved 
by extracting heat through parametric damping or ac- 
tive feedback [l|, Theoretical calculations suggest 
that oscillators with a large thermal occupation number 
(fcsT ^ hijjm) can be cooled to be close to their ground 
state, if they have high enough quality factors Once 
the ground state is approached, many interesting stud- 
ies of macroscopic quantum mechanics can be performed, 
e.g. teleporting a quantum state onto mechanical degrees 
of freedom , creating quantum entanglement between a 
cavity mode and an oscillator Q and between two macro- 
scopic test masses Most proposals involve the oscil- 
lator position linearly coupled to photons, in which case 
the quantum features of the oscillator, to a great extent, 
are attributable to the quantization of photons. In or- 
der to probe the intrinsic quantum nature of an oscilla- 
tor, one of the most transparent approaches is to directly 
measure its energy quantization, and quantum jumps be- 
tween discreet energy eigenstates. Since linear couplings 
alone will not project an oscillator onto its energy eigen- 
states, nonlinearities are generally required 0, SQ- For 
cavity-assisted optomechanical systems, one experimen- 
tal schemejproposed in the pioneering work of Thomp- 
son et al. []|, is to place a dielectric membrane inside a 
high-finesse Fabry-Perot cavity, forming a pair of coupled 
cavities [l^l ■ If the membrane is appropriately located, a 
dispersive coupling between the membrane position and 
the optical field is predominantly quadratic, allowing the 
detection of mechanical energy quantization. 

In this letter, we show that in the experimental setup 
of Thompson et al., the optical field also couples linearly 
to the membrane. Due to finiteness of cavity finesse (ei- 
ther intentional for readout or due to optical losses), this 
linear coupling introduces quantum back-action. Inter- 
estingly, it sets forth a simple standard quantum limit, 
which dictates that only those systems whose cavity- 




FIG. 1: The left panel presents the schematic configura- 
tion of coupled cavities in the proposed experiment [J] . The 
right panel shows optical modes and we denote refiectivity and 
transimisivity of the optical elements by and ti {i = 1, 2, m). 



mode decay rates are smaller than the optomechanical 
coupling rate can successfully resolve energy levels. We 
will further show that a similar constraint applies univer- 
sally to all experiments that attempt to probe mechanical 
energy quantization via parametric coupling with exter- 
nal degrees of freedom (either optical or electrical). 

Coupled Cavities. — Optical configuration of coupled 
cavities is shown in Fig.[TJ Given the specification in Ref. 
[l[ , transmissivities of the membrane and end mirrors are 
quite low, and thus a two-mode description is appropriate 
[lol . ITU , with the corresponding Hamiltonian 

ii = hujjnif +P^)/2 + huJo{a''a + b^b) - hu;s{a^b-\-b^a) 
+ hGoq{a^a-Pb)+n^^t+ne- (1) 

Here q,p are normalized position and momentum of 
the membrane; a, b are annihilation operators of cav- 
ity modes in the individual cavities (both resonate at 
loq); Us = tmc/L is the optical coupling constant for 
a and b, through transmission of the membrane 
Go = 2\/2ijJoXq/L is the optomechanical coupling con- 
stant with L denoting the cavity length and zero-point 

Hcxt and correspond to 



motion Xq = ^Jh/ {2muom)\ 
the coupling of the system to the environment and quan- 
tify the fluctuation and dissipation mechanism. By intro- 
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ducing optical normal modes, namely the common mode 
c = (a 4- b)/V2 and differential mode d= {a — 6)/\/2, 

+ t{y^cc''cin + y^d^di^-}i.c.)+n^/h (2) 

where w± = ujoiajs and in the Markovian approximation 
H.cxt is written out explicitly in the second line (with jc,d 
denoting decay rates and H.c. for Hermitian conjugate). 

Before analyzing the detailed dynamics, here we follow 
Thompson et al. jl| and Bhattacharya and Meystre 
by assuming ujm <?C oJs and Go ^ — lu-\ = 2ujs, 
analogous to the dispersive regime in the photon-number 
counting experiment with a superconducting qubit [l^, 
13|. This allows us to treat hGoq{c'd + rftc) as a per- 
turbation and diagonalize the Hamiltonian formally. Up 
to Gq/{2lUs)'^ , the optical and optomechanical coupling 
parts of the original Hamiltonian can be written as 



n/h: 



2u, 



(3) 



At first sight, frequency shift of the eigenmodes 6 and e 
is proportional to (f' . Since frequency separation of two 
normal modes is 2 cj^ 3> 7c, d, they can be independently 
driven and detected. Besides, with jc.d < w™, only av- 
eraged membrane motion is registered and = -|- 1/2 
with N denoting the number of quanta. Therefore, previ- 
ous authors had concluded that such a purely dispersive 
coupling allows quantum non-demolition (QND) mea- 
surements of the mechanical quanta. 

However, the new eigenmodes 6 and e are given by 

o = c - [{God)/{2uj,)]q, e^d+ [{Gqc)/ {2io,)]q. (4) 

If we pump c with classical amplitude c and left d in vac- 
uum state, the detected mode 6 will have a negligible lin- 
ear response. However, the idle mode e w [Goc/{2ujs)]q, 
which is dominated by linear coupling. If we choose to 
drive d, the role of 6 and e will simply swap. Such lin- 
ear coupling can potentially demolish the energy eigen- 
states that we wish to probe. We can make an order-of- 
magnitude estimate. The optomechanical coupling term 
in Eq. at the linear order, reads Goq{cd + c*d''). Ac- 
cording to the Fermi's golden rule, it causes decoherence 
of energy eigenstate near the ground level at a rate of 



'dec 



where we have assumed that c is on resonance, and 



S^^ = jdte'^\dlt)d){0)) = 2^a/[{Lo - 2lo, 



(5) 



(6) 



On the other hand, from Eq. ([3]) and linear response the- 
ory IJ], the measurement time scale to resolve the en- 
ergy eigenstate (i.e. measuring N with a unit error) with 
a shot-noise limited sensitivity is approximately given by 

r„, « [^lu:l/{Gt,\-cn~Sm = 2ul^J{Gl\-c\% (7) 



where S'c(O) is the spectral density of c at zero frequency. 
Requiring t,„ < Tdcc yields 



ilcldlGl) < 1. 



(8) 



In the case when transmissivity of end mirrors ti ^ t2 = 
to, we have 7c = 7d = ct^l{2L). Defining the cavity 
finesse as J- = T^/t^, the above inequality reduces to 
\/{J-Xq) < 8a/2. Therefore, to probe mechanical energy 
quantization, it requires a strong-coupling regime (c.f. 
Eq. ([5])), or equivalently, for such an optomechanical sys- 
tem, zero-point mechanical motion Xq to be comparable 
to linear dynamical range \/ T oi the cavity. 

We now carry out a detailed analysis of the dynamics 
according to the standard input-output formalism [l5| . 
In the rotating frame at the laser frequency lj^ , the non- 
linear quantum Langevin equations are given by 

q = ujyyi p, (9) 

P = -^mq~ lrnP~GQ{c^d + d)c)+ith, (10) 

c = -7cC- iGoqd-h y27^Cin, (11) 

d = -{'^d + 2ii^s) d~ iGoqc+ y/2ri2d{n- (12) 

Here the mechanical damping and associated Brownian 
thermal force ^th origin from Ti.^ under the Markovian 
approximation. These equations can be solved perturba- 
tively by decomposing every Heisenberg operator d into 
different orders such that a = a + e d^^^ -|- e^d^^^ -I- ©[e^]. 
We treat G{)/{2ujs), vacuum fluctuations ^/27c cl^ and 

^j2qdS^ (simply denoted by ^flrfchn and ^J2qdd\n in 
later discussions) as being of the order of e (e ^ 1). 

To the zeroth order, c = y^2Io/{"fcf'-ujo) with Iq denot- 
ing the input optical power and d = 0. Up to the first 
order, the radiation pressure term reads Goc[d'-^^ + d'^^^^ 
(c is set to be real by choosing an appropriate phase ref- 
erence). In the frequency domain, it can be written as 

_ 2y/^ Go c[(7d - iuj)vi - 2ujsV2] + ^Glc^Usg 
'•^ (oi + 2^^. + «7d)(w - 2ws + i7d) 

where vi , V2 and q are Fourier transformations of i'l (t) = 
{din + dl)/V2, V2{t) = {din-dl)/{iV2) and q{t) respec- 
tively. The part, containing vacuum fluctuations, is the 
back-action Fba, which induces the quantum limit. The 
other part proportional to q is the optical-spring effect. 
Within the time scale for measuring energy quantization, 
of the order of 7~^ (<g; 7^^), the positive damping can 
be neglected but the negative rigidity has an interesting 
consequence — it modifies tUm to an effective uJeS (< ^^m)- 
Correspondingly, position of the high-Q membrane is 

q{t) = <irn + A^Jodt' SmUJ,ff{t~t')[FBA{t')+^th{t')] (14) 



with A y^uJm/uJcS- The free quantum oscillation 
(jm ~ A ((^0 coswoff t + po sinweff t) and go and pQ are the 
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initial position and momentum normalized with respect 
to ^Jh/{muJcs) and y^hmcueg. 

The dispersive response is given by the second-order 
perturbation 0(6^] . Adiabatically eliminating rapidly os- 
cillating components and assuming ^ which can 
be shown to maximize the signal-to-noise ratio, we obtain 



(15) 



Here Gcg = A Gq and N{t) = Nq + AN{t) contains 
the number of mechanical quanta A^o = (^Zo +Po)/2 ^-nd 
the noise term AN{t) due to the back-action and ther- 
mal noise. To read out N{t), we integrate output phase 
quadrature for a duration r. According to the input- 
output relation Cout + Cjn = \/27c c, the estimator reads 

f(r) = J,'dt[Mt) ~ ciV(t)/(V^^,)], (16) 

where U2 = (cin — cj^J/ {i\/2). For Gaussian and Marko- 
vian process, the correlation function {c2{t) c\{t')) = 
S{t — t')/2. For typical experiments, the thermal occu- 
pation number fith = ksT / {huj„i) is much larger than 
unity, and (Cth(0 fth(i')) ~ 2^rnnth5{t - t'). Through 
evaluating the four-point correlation function of back- 
action noise and ^tii(i) in {AN{t)AN{t')), wc obtain the 
resolution AA^ as a function of r 



AA^^ 



5 / 7rfGgffcV 

6 V 2V2tj2 



5 f -frnksTT 

6 



V2h 



(17) 

In order to successfully observe energy quantization, 
the following conditions are simultaneously required: (i) 
the resolution AA^^ should have a minimum equal or less 
than unity, (ii) this minimum should be reachable within 
r that is longer than the cavity storage time 1 /jc (which 
in turn must be longer than the oscillation period 1 /wcff 
of the membrane), (iii) the system dynamics should be 
stable when taking into account optical rigidity which is 
approximately equal to GqC^/uJs for ujm Ws- 

Specifically, the standard quantum limit in condition 
(i), set by the first two terms in AA^^, gives jdd/Glff < 
1, or equivalently {jcJd/Gl) < A^. If we neglect the op- 
tical spring effect (A = 1), wc simply recover Eq. A 
strong negative optical rigidity (ujcS ^ ^m, i.e. A ^ 1) 
can significantly enhance the effective coupling strength 
and ease the requirements on optomechanical properties. 
However, a small UeS also makes the system susceptible 
to the thermal noise. Taking account of all the above con- 
ditions, the optimal ujcff = (^m\/nth/Qm with mechani- 
cal quality factor Qm = i^m/jm, and there is a nontriv- 
ial constraint on the thermal occupation number, which 
reads (fith/Qm) < [Gg/(tJ^7c)]^/^. 

For numerical estimate, we use a similar specification 
as given in Ref. [H but assume a slightly higher me- 
chanical quality factor Qm, lower environmental temper- 
ature T and lower input optical power Iq such that all 




FIG. 2: The resolution AA^ for measuring mechanical energy 
quantization depending on the integration duration r with 
total noise (Solid) and quantum noise only (Dashed). 



mentioned conditions are satisfied. The parameters are 
the following: m = 50 pg, cj„i/(27r) = 10^ Hz, Q,n = 
3.2 X 10^, A = 532 nm, L = 3cm, r„ = 0.9999, T = 
6 X 10^ T = 0.1 K and lo ~ 5 nW. The resulting resolu- 
tion AN is shown in Fig. [21 and we are able to resolve 
single mechanical quantum when r w 0.1 ms. 

Even thought we have been focusing on the double- 
sided setup where ti « t2, the quantum limit also exists 
in the single-sided case originally proposed in Ref. [H. 
Ideally, a single-sided setup consists of a totally reflected 
end mirror and the vacuum fluctuations only enter from 
the front mirror. Therefore, the quantum noises inside 
two sub-cavities have the same origin but different opti- 
cal path. Through similar input-output calculations, we 
find that if laser detuning is equal to iu;^, the quantum 
noises destructively interfere with each other at low fre- 
quencies, due to the same mechanism studied in great 
details in Ref. , achieving an ideal QND measure- 
ment. However, in reality, the end mirror always has 
some finite transmission or optical loss which introduces 
uncorrelated vacuum fluctuations. As it turns out, the 
quantum limit is similar to Eq. only with 7c, d re- 

placed by the damping rate of two sub-cavities. 

General Systems. — Actually, the standard quantum 
limit obtained above applies to all schemes that attempt 
to probe mechanical energy quantization via parametric 
coupling. Let us consider n mechanical modes parametri- 
cally coupled with n' normal external modes, describable 
by the following Hamiltonian 



(18) 



n' n 



Here Greek indices identify mechanical modes and Latin 
indices identify external modes; flu and tui are eigenfre- 
quencies; QutPu are normalized positions and momenta; 
di are annihilation operators of the external degrees of 
freedom; Xijv = Xjii' are coupling constants. Similarly, 
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we focus on the regime where |xyv| ^ 



(disper- 



sive) and ^ \uJi 



(adiabatic), and obtain 



7i = ^;jf7,(<72+p2)/2+^ft^^'al6,+?i:ext+H?, (i9) 



where, up to X%u/W^ ~ ^j?^ 



... LOi-UJj 



(20) 



In order to have quadratic coupUngs between a pair of 
external and mechanical modes, b\ and q\ for instance, 
we require that Xwv = and xuu = Xui^iu^ ^-iid then 



91- 



1^1 



(21) 



However, there still are linear couplings which originate 
from idle modes. This is because, up to Xijiy/\^i ~ 



Xi]iaj 



qi [i ^ 1). (22) 



where ai is replaced with its classical amplitude di, for 
ai » fli. From Eq. ([2T|) and ([22]) . both linear and dis- 
persive couplings are inversely proportional to — wij. 
Therefore, we only need to consider a tripartite system 
formed by qi, 6i and 62 which is the closest to 61 in 
frequency. The resulting Hamiltonian is identical to Eq. 
(PI), and thus the same standard quantum limit applies. 

Conclusion. — We have demonstrated the existence of 
quantum limit for probing mechanical energy quantiza- 
tion in general systems where mechanical modes para- 
mctrically interact with optical or electrical degrees of 
freedom. This work will shed light on choosing the ap- 
propriate parameters for experimental realizations. 
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